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ABSTRACT 


The  hydromagneti'c  model  of  the  evolution  of  a  plasma  puff  in  a 
field  described  by  Waelbroeck  has  been  used  to  analyze  the  character¬ 
istics  of  the  plasma  produced  by  the  APL  conical  pinch  gun.  The  analy¬ 
sis  is  based  on  fitting  the  signals  induced  in  single  turn  magnetic  pick¬ 
up  loops  by  the  passing  plasma.  In  the  cases  studied  good  fits  to  the 
signals  have  been  achieved  and  the  characteristic  parameters  of  the 
plasma  found.  This  analysis  applied  to  plasma  generated  by  the  APL 
theta-pinch  gun  has  not  been  successful.  In  the  two  cases  studied,  the 
lack  of  good  fits  is  attributed  to  the  presence  of  a  sizeable  amount  of 
trapped  field  in  the  puff  which  is  present  along  the  entire  flight  path 

through  the  B  field, 
z 


The  John*  Hopkint  Univeriily 
APPLIED  PHYSICS  LABORATORY 
Silvtf  Spring,  Mtrylond 


TABLE  OF  CONTENTS 


I.  Introduction  . 

II.  Theory  . 

III.  Magnetic  Loop  Signals  Induced  by  the  Plasma 

IV.  Determination  of  Parameters  from  Magnetic 

Signals  . 

V.  Analysis  . 

VI.  Computer  Analysis  . 

VII.  Results  . 

VIII.  Discussion  . 


1 

3 

16 

21 

26 

29 

41 

43 


Appendices 

I.  Solution  of  Equation  1  Ob 

II.  Solution  of  Equation  1 2 

III.  Solution  of  Equation  22 

IV.  The  Compensated  Magnetic  Loops 

V.  Computer  Listings 

References 


44 

46 

47 
50 
52 

62 


v 


Th*  Johm  Hopkini  Univanity 

APPIICO  PHYSICS  LABORATORY 

Sllvof  Spring,  Maryland 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Experimental  Arrangement  .  .  .  12 

2  Magnetic  Loop  and  Field  Geometry  .  .  17 

3  Domain  of  Allowable  Values  of  (f*,  A514)  •  •  24 

4  Typical  Magnetic  Loop  Signals  .  .  .  33 

5  Comparison  of  Magnetic  Loop  Signals  with 

Analysis  (points),  Condition  I  .  .  35 

6  Comparison  of  Magnetic  Loop  Signals  with 

Analysis  (points),  Condition  II  .  38 


Th«  Johm  Hopkim  Univ«f»i!y 

APPLICO  PHVtlCft  IAAOAATOAY 

Silver  Spring,  M*ryUnd 


I.  INTRODUCTION 

This  report  is  an  elaboration  of  the  work  conducted  by  the  fusion  re¬ 
search  group  (Ref.  1)  at  Fontenay-aux-Roses  during  the  period  1960-1964 
concerning  the  analysis  of  plasma  flow  through  a  longitudinal  magnetic 
field.  This  report  includes  a  more  comprehensive  discussion  of  the  model 
which  was  presented  by  the  French  scientists  as  well  as  the  analysis  of 
plasma  flow  made  by  the  Plasma  Dynamics  Group  at  the  Applied  Physics 
Laboratory.  The  computer  analyst  developed  at  the  Laboratory,  is 
also  discussed  in  some  detail. 

The  model  describes  the  plasma  by  means  of  hydromagnetic  equations 
with  some  simplifying  assumptions.  It  is  shown  that  the  model  is  able  to 
extract  the  gross  features  of  the  plasma  (such  as  temperature,  density, 
center  of  mass  and  local  velocity)  from  magnetic  loop  signals  induced  by 
the  plasma  while  passing  along  the  field.  The  hydromagnetic  description 
makes  it  immediately  obvious  that  this  analysis  is  applicable  to  a  rather 
restricted  class  of  plasmoids.  Nevertheless,  in  the  original  work  (Ref. 

1)  it  was  found  that  whenever  the  model  could  be  made  to  fit  the  experi¬ 
mental  data,  the  theoretical  model  predictions  and  those  of  more  direct 
diagnostic  techniques  were  in  excellent  agreement.  Therefore,  this  theory 
is  considered  very  useful,  albeit  to  a  restricted  family  of  plasmas. 

Both  the  French  workers  and  those  at  APL  analyzed  plasmas  which 
were  produced  by  a  conical  pinch  gun.  The  gun  coils  used  in  both  cases 
were,  in  most  respects,  identical.  However,  whereas  the  French  gun  was 
powered  by  a  1 0  ^xf  capacitor  bank,  the  APL  gun  used  three  14-/if  capaci¬ 
tors  and  had  a  nominal  frequency  of  80  kc.  The  maximum  allowable 
charging  voltage  was  limited  to  14-kv  by  the  current  capacity  of  the  igni- 
tron  switches. 


-  1  - 


JWWTH  . 

“If® 


Th«  John i  Hopkint  Unlvartlly 
APPLIED  PHVtICt  LABORATORY 

Sllv*r  Spring,  M»ryl«nd 


It  is  hoped  that  this  analysis  can  be  extended  to  the  diagnosis  of  a 
plasma  puff  in  a  multipole  magnetic  field.  The  work  discussed  here  and  the 
investigation  of  this  more  complicated  magnetic  configuration  are  being  sup¬ 
ported  by  the  National  Aeronautics  and  Space  Administration,  Physics  of 
Fluids  Branch. 
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II.  THEORY 


The  evolution  of  a  plasma  puff  as  it  passes  through  a  longitudinal  mag¬ 
netic  field  has  been  described  by  Waelbroeck  et  al  in  Ref.  1.  The  hydro- 
magnetic  equations  have  been  used  with  the  following  assumptions: 

1.  The  longitudinal  guide  field  does  not  impede  axial  motion  of  the  plasma, 
but  it  does  eliminate  radial  diffusion.  The  density  is  sufficiently  low  that 
diffusion  across  the  magnetic  field  is  not  appreciable  during  the  time  of  in¬ 
terest. 

2.  At  any  time,  t,  and  for  a  given  axial  position,  the  plasma  characteris¬ 
tics  are  independent  of  the  radial  coordinate,  r,  up  to  the  boundary  of  the 

plasma,  r  . 

P 

3.  The  collision  frequency  is  sufficiently  high  that  the  perpendicular  and 
parallel  temperatures  (with  respect  to  the  direction  of  the  guide  field)  are 
equal  at  every  point,  i.  e.  ,  Te|j  =  T^,  T.  |j  =  T. 

4.  The  plasma  is  electrically  neutral  (n  .  £  =  n^)  and  each  fluid  element, 
either  electron  and  ion  gas,  has  the  same  average  velocity  (V  =  V.),  where 

C  1 

n^  ,  n^  and  £  are  the  ion  density,  electron  density  and  ionic  charge,  respec¬ 
tively.  This  assumption  eliminates  any  longitudinal  current. 

5.  The  magnetic  pressure  and  plasma  pressure  are  in  equilibrium  at 
every  instant. 

6.  At  every  .instant,  the  linear  density  of  the  ions  in  the  plasma  with 
respect  to  the  center  of  mass  of  the  puff  is  Gaussian,  i.  e.  , 


N. 

i 

V^T  X(t) 


exp 


(1) 


where  N.  is  the  total  number  of  ions  in  the  puff,  z  is  the  coordinate  rela¬ 
tive  to  the  center  of  mass  of  the  puff  and  A(t)  is  the  characteristic  width 
of  the  Gaussian  distribution. 
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The  momentum  nsfer 
can  be  written  in  the  .  uiowing 


equation  with  scalar  pressure  (see  Ref. 
form: 


2) 


p  ^7  -  VV  •  (p  V)  +  V  (n.m.V.V.)  +  V  (n  m  V  V  ) 
r  9  t  ^  1111  eeee 


g  E  +  (j  x  B)  -  V  p  -  pVcp  , 


(2) 


where  V.  and  V 

l 

of  mass. 


e 


the  ion  and  electron  velocities  with  respect  to  the  center 


m.  and  m 
l  e 


are  the  corresponding  masses  of  these  particles. 


cr  =  (nr  -  n  )  e 
l  e 

(charge  density) 

(3a) 

p  =  n.m.  +  n  m 
r  ii  e  e 

(mass  density) 

(3b) 

V  =  —  (n.m.  V.  +  n  m  V  ) 
p  ill  e  e  e 

(macroscopic  velocity) 

(3c) 

3  =  (nr  V.  -  n  V  )  e 

1  P 

(current  density) 

(3d) 

i  i  e  e 


p  =  the  total  scalar  pressure  (electron  plus  ion  pressure) 
B  =  the  applied  magnetic  field,  and 
cp  =  the  gravitational  potential. 


Neglecting  the  gravitational  term  and  the  electron  mass  compared  to 
the  ion  mass,  Eq.  2  becomes: 


n.m. 
i  i 


-  V.  V  •  (n.m.  V.)  +  V  • 

l  ill 


(n.m.V.V.)  = 
1111 


a  E  +  (J  x  B)  -  Vp  . 


(4) 


But, 


V(n.m.  V.V. ) 
1111 


n.m.V.  •  v  V.  +  V.V(n.m.V)  . 

ill  i  i  ill 


(5) 
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It  is  interesting  to  note  that  this  equation  describes  the  flow  of  a  simple 

fluid  with  the  thermal  energy  density,  3/2  n.k  (T  +  T.),  under  the  force 
1 

-  —  v  p  (where  k  is  the  Boltzmann  constant), 
n. 

1 

The  basic  parameters  are  defined  as  follows: 

1.  Xq  and  X(t)  are  the  characteristic  widths  of  the  puff  along  the  axis 
at  t  =  o  and  at  any  time  t,  respectively.  It  should  be  noted  that  X(t)  is  a 
function  of  time  only.  Time  zero  is  defined  as  the  time  at  which  the  puff 
begins  to  expand  axially  and  is  not  necessarily  the  time  of  ejection  from  the 
gun. 

2.  vq  is  the  velocity  of  the  center  of  mass  of  the  puff.  v(Z,  t)  is  the 
velocity  of  any  element  of  the  fluid  at  a  distance,  Z  centimeters,  from  the 
initial  position  of  the  center  of  mass  and  at  a  time,  t  seconds,  after  the 
expansion  begins. 

3.  Bq  is  the  guide-field  intensity. 

4.  r^  is  the  radius  of  the  flux  conserver  on  which  the  solenoid  is  wound. 

5.  r  is  the  radius  of  the  flux  tube  in  which  the  plasma  travels  and  will 

P 

vary  with  position  and  time,  r  is  the  radius  of  this  flux  tube  when  the 
plasma  is  not  present.  The  total  magnetic  flux  inside  the  flux  conserver 
is  constant. 
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The  volume  density  of  the  ions,  expressed  in  a  form  similar  to  Eq.  1, 
is  given  by 


N, 


n. 


Jn  X(t) 


exp  - 


r-  -n  2 

z 


(8) 


where  N  is  the  number  of  ions  in  an  infinitely  long  cylinder  of  unit  cross- 
sectional  area.  Using  this  expression  in  the  equation  of  continuity, 


3  n. 

-T7 -  +  “  (n.V.)  =  0  , 

3t  O.z  1  1 


(9) 


one  obtains 


3  n. 
_ i 

3t 


n. 

l 


_£Z_ 


X 


I 

X  ; 


dX 

dt 


(10a) 


and 


n.V.  =  -  X 

l  l 


2z2n. 


n. 

l 


L  Xc 


d.z 


(10b) 


It  can  be  shown  (see  Appendix  I)  that  the  solution  for  this  velocity  relative 
to  the  center  of  mass  which  is  anti- symmetric  for  1  z  is 


V. 


z  X 


(ID 


The  equation  of  motion,  Eq.  7,  now  becomes 


n.m.zX 
3p  ii 


3z 


+ 


X 


=  0 , 


(12) 
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which  has  the  solution  (see  Appendix  .11)  expressed  by 


p  =  —  n.m.  X  X" 
2  1  1 


Since  it  has  been  assumed  that  T||  =  Tj^  for  each  component  gas, 
the  ideal  gas  law  is  applicable,  i.  e.  , 


p  =  n.k  T.  +  n  k  T  =  n.k  (T.  +  T  )  , 
liee  lie 


and  thus 


k  (T  +  T.)  =  -  m.  XX" 
e  l  2  l 


From  this  it  may  be  seen  that  the  temperature  is  a  function  of  time  only  and 
therefore  is  uniform  in  space.  This  implies,  for  a  Gaussian  density  distri¬ 
bution  and  a  plasma- field  equilibrium  at  each  instant,  that  the  plasma  radius, 

r  ,  varies  with  z. 

P 

The  total  energy  of  the  plasma  puff  will  be  considered  to  exclude  the 
effect  of  external  forces.  This  total  energy  is  composed  of  the  thermal  en¬ 
ergy  of  the  ions  and  the  electrons  plus  the  kinetic  energy  of  the  ions'  axial 
motion.  The  kinetic  energy  of  the  axial  motion  of  electrons  is  negligible 
since  and  rm  »  m^  .  The  energy,  W,  is  defined  as 


W  =  |  f”  f  P  n.k  (T  +  T.)  2  nr  drdz  + 
2  J-m  i  e  i 


1  00  n  P  P 

—  f  f  n.m.V.1 2  2  ij  r  drdz 

2  J-cc  Jo  ill 
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The  ideal  solution  would  be  to  integrate  Eq.  16  over  the  radius  of  the  plasma 
and  thus  calculate  the  energy  per  unit  area.  Equating  the  pressure,  Eq.  14, 
to  the  energy  density  would  then  give  a  useful  relationship.  Unfortunately,  r 
is  not  known  as  a  function  of  z  and  thus  the  integral  cannot  be  evaluated 
exactly.  However,  if  only  an  undefined  low  ft  plasma  is  considered,  the 
plasma  radius  can  be  said  to  vary  slowly  with  z.  Equation  16  may  now  be 
written  as 


W 


n,m.  V.2 

ill 


dz 


(17) 


Combining  Eqs.  11,  15,  and  17  and  integrating  the  result  yields 


_W_  =  JL  m.  N.  XX"  +  1  m.  N  X1  ,  (18) 

7rr  2  4  i  A  4  l  A 

P 


Within  the  limits  of  the  restricted  j3,  this  expression  is  equal  to  a 
constant  which  can  be  evaluated  at  time  zero,  i.  e.  , 


and  hence 


at  t  =  o,  X  =  X  and  X'  =  o  , 
o 


-  NA  k  (T  +  T.)  _  . 
2  A  e  i  o  * 


The  approximation  of  low  /3  must  be  examined  more  carefully.  The 
J3  of  the  plasma  is  defined  in  the  usual  way  as 


0=  p(Be2/2pQ) 


(19) 
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where  B  is  the  field  inside  the  plasma  boundary  and  B'  is  the  field  on  the 
vacuum  side  of  this  boundary.  The  external  field,  B',  is  approximately 
equal  to  the  unperturbed  field,  B  ,  as  long  as  the  plasma  radius  is  small 
compared  with  the  radius  of  the  solenoid  (this  is  usually  the  case).  Then, 
from  the  conservation  of  flux, 

it r2  B  =  it r2  B  , 
p  o  o 


and  Eq.  19  takes  the  form 


£ 


or 


—  =  (  1  -  0)*  • 
r 

P 


(20) 


When  (3  =  0.  6,  the  ratio  of  the  flux  tube  radius  to  the  plasma  radius  is  about 
0.  8.  Thus  even  for  relatively  high  j3  plasma,  the  plasma  radius  changes 
slowly  with  z  and  the  approximation  used  in  Eq.  17  is  a  fairly  good  one. 

The  following  reduced  variables  are  now  defined: 

A  =  4^  ,  (21a) 

Ao 

and 

t 

r  =  T  • 


where 


a  = 


m.X2 
_ 1  o 

6k  (Te  +  Ti)0 


(21b) 
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Differentiating  Eq.  21a  yields 


_  d_X  _  dA  dT _ o  dA 

dt  Ao  dr  dt  "  £  d r  ' 


da  X  _  ^o_  _£A 

.  o  2 


£  dr 


Equation  18  can  now  be  written  in  terms  of  these  variables  as 


\  a  x  r  a 

3  .  o  daA  4.  1  „  o  dA  _  3  XT  /rT,  ,  m  % 

-  N .  m.  A  —5-  - -  +  7  N  m.  ——  — —  -  —  N.k  (T  +  T.) 

4  A  1  dT=  4  A  1  d  T  2  A  e  1  o 


which  reduces  to 


3A^4  ♦ 


This  differential  equation  expresses  the  conservation  of  energy.  The  solu¬ 
tion  of  this  non-linear  differential  equation,  as  given  in  Appendix  III,  is 


T  =  (A*  +  2)  <A^  -  1)*, 


=  (1  -  A  *)  . 


These  equations  combine  with  Eq.  15  to  yield 


k  (T  +  T.)  =  3  k  (T  +  T.)  A 


e  1 


e  1  o  d  t‘ 


and  thus 


(T  +  T.) 
e  1 

(T  +  T.) 


rx  1 
_ 0_ 

X  (t) 

-  — 
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The  information  which  has  been  gathered  thus  far  about  the  plasma  is 
not  in  a  useful  form.  These  equations  must  be  expressed  in  the  laboratory 
coordinate  system  rather  than  in  the  center-of-mass  system  and  they  must 
be  related  to  experimentally  measurable  quantities. 

Beginning  with  the  coordinate  transformation,  it  can  be  seen  that  the 
laboratory  axial  position  (measured  from  the  center  of  the  gun  coil  to  any 
location,  z,  relative  to  the  center  of  mass  of  the  puff,  see  Fig.  1)  is 
given  by 

z  =  z  +  v0tQ  +  VQt  =  z  +  ZQ  +  vQt  ,  (26) 

where  ZQ  =  vQt0  is  the  laboratory  position  at  which  the  puff  begins  to  expand. 
Thus. 

z  =  Z-  Z  -vt  •  (27) 

o  o 

The  fact  that  expansion  does  not  take  place  immediately  can  probably 
be  accounted  for  as  follows: 

1.  Ions  formed  downstream  are  ejected  earlier  than  other  ions  but  with 
a  lower  velocity.  Thus,  the  ions  tend  to  bunch  initially. 

2.  Since  the  model  does^not  include  the  possibility  of  a  trapped  field  which 
inhibits  expansions,  the  analysis  is  applicable  only  after  this  field  has  been 
dissipated. 


Some  of  the  expressions  derived  earlier  may  now  be  rewritten  as 


X(t)  X(t)  A  A 


SOLENOID  WINDING 


Fig.  1  Experimental  Arrangement 
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Defining 


and 


:v 


o 


a 


X 


o 


2 

m .  v 
1  © 


6k(T  +T.) 
e  1 


o-J 


Z  -  Z 


i  - 


o 


X 


o 


Eqs.  1  and  8  become,  respectively, 


and 


exp 


2 

Z '  -  aT  ~ 

A 

2 

Z'  -  aT 

A 


The  transverse  energy  per  unit  length  of  the  plasma, 


W 


is 


N  k  (T  +  T.)  = 
l  e  i 


N,  k  <T  +  T.) 
i  e  i 

‘Ai  A  Q  A 


exp 


“1  2 


Z'  -  aT 


But  from  Eq.  25  we  have 


and  thus 


where 


(T  +  T.) 
e  i 


(T  +  T.) 

e  i  o 

- — ? — 

A  5 


N.  m.  v 

w  =  -1 — I — - —  f  (A.,  o>  Z1)  , 

-L  6  \/ji  a  £ 


f  (A,  a,Z) 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 
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(3  may  now  be  expressed  in  terms  of  the  linear  transverse  energy 
density.  Since 


2mdp 


2  77  r  jj  p 
p  o 

77  r2  B2 
p  o 


P  77  r  =  W  .  , 
p  Jl± 


then 


”  77r2  B2 

P  o 


=  2A  | 


where  A  is  defined  as 


H  W,  , 

A  -  12 _ LL 

A  '  77r2  B2 
o  o 


u  N.m.v 
oiio 

6  77  3^2  a  £  r2  B2 
^  o  o 


f  (A,o  ,Z  ) 


But  (3  has  been  defined  by 


---i-  i- 


From  this  it  follows  immediately  that 


=  Ja*  +i  -  a  . 

B  r  v 

o  L  pJ 


14  - 


Ih»  John*  Hoptiln*  Unlv«f»ity 
APPLIED  PHVSIOS  LABORATOKV 
Silvor  Spring.  M*fyl*nd 


From  Eq.  (32)  the  ion  density  may  now  be  written  as 


N 


N. 


1 


n.(Z',r)  =  —  =  — 

1  jrr2  W3X  r»A 


exp  - 


Z'-or 


L  A  J 


or 


Ni 

n (Z, t )  =  - 

173  /3  X  r3 
o  o 


o  p 

(v^a+l  -A)  exp  - 


r  i  -i  a 
Z  -QiT 

-A-a._ 


(41) 


The  local  plasma  velocity,  v(Z  ,r)  may  also  be  obtained.  It  is  the 

sum  of  the  center-of-mass  velocity  and  the  ion  velocity  (with  respect  to  the 

*  z\  * 

center  of  mass),  i.  e.  ,  v(Z  ,r)  =  V  +  V.  =  V  +  — —  ,  which  becomes 

o  i  o  X 


or 


,  X 
v(Z.T)  =  f 

v(z'.r)  *  j 


(42a) 


(42b) 


It  can  be  seen  from  the  above  discussion  that  the  evolution  of  the  plas¬ 
ma  puff  can  be  described  by  six  parameters:  v  ,  (T  +T.)  ,  Z  ,  X  ,  N.  and 
.  .  oeioooi 

r^  or,  equivalently,  by  a,  £,  Zq,  Xq,  and  rQ.  The  remainder  of  this 

report  will  discuss  means  of  finding  these  parameters  from  an  analysis  of 
magnetic  loop  signals. 
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III.  MAGNETIC  LOOP  SIGNALS  INDUCED  BY  THE  PLASMA 

The  work  of  Waelbroeck  et  al  (Ref.  1)  differs  slightly  from  the  APL 
work  in  that  the  French  investigators  used  a  special  type  of  magnetic  probe, 
called  a  compensated  magnetic  loop,  while  simple  single-turn  pick-up  loops 
have  been  used  at  APL.  The  compensated  loops  are  discussed  in  Appendix 
IV  and,  as  can  be  seen,  have  some  definite  advantages  over  the  system  used 
here.  However,  since  one  of  the  primary  purposes  of  the  present  work  is 
to  extend  this  analysis  to  the  multipole  magnetic  configuration  where  only 
simple  loops  can  be  used,  the  convenience  of  the  compensated  loops  is 
sacrificed. 

Figure  2  is  a  schematic  representation  of  the  drift  field  and  a  simple 
magnetic  loop.  The  period  of  the  field  is  very  long  compared  to  the  transit 
time  of  the  plasma,  and  the  field  is  spatially  uniform  inside  the  flux  con- 
server.  The  plasma,  effectively  excludes  part  of  the  field  from  the  volume 
which  it  otherwise  occupies.  In  the  presence  of  the  flux  conserver,  this  ex¬ 
cluded  flux  is  distributed  over  the  volume  between  r  and  the  conserver 

P 

radius,  r  .  For  a  plasma  whose  length  is  long  compared  to  r  ,  this  dis- 

O  O 

tri'oution  can  be  (considered  uniform. 

Assume  Bq  is  the  drift  field  strength  in  the  absence  of  plasma, 

-Bp  is  the  field  generated  by  currents  internal  to  the  plasma, 
and 

is  the  increase  in  the  field  strength  external  to  the  plasma 
caused  by  the  presence  of  the  plasma. 

With  the  previous  assumption  that  the  plasma  puff  contains  no  trapped  field, 
the  increase  in  flux  between  r^  and  r^  may  be  equated  to  the  decrease  in 
flux  inside  the  plasma.  If  Ap  represents  this  flux  change,  then 
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A<pp  =  2ir  J  Bp  r  d  r  =  2  ff  f  Berdr. 

O  Tp 


(43) 


Since  Be  is  assumed  to  be  uniform  between  Tp  and  rg  ,  Eq.  43  becomes 


A<n. 


B 


6  ir(  rl-rp 


(44) 


The  flux,  <p(t),  at  any  time  t  encircled  by  our  magnetic  loop  (shown 
in  Fig.  2)  may  be  expressed  as 


<P<t>  =  irv3  (Bc-B  )  +  ir(r?  -  r?)  (B0  +  Be> 

I _ P  Pl  I  P _ I 

flux  inside  flux  outside  plasma 

plasma  encircled  by  magnetic 

loop 


=  -  A<pp  +  ff(i*®  -  rp  Be  + 


B0 
nvt  B0 


The  emf  generated  is  then  given  by 


€ 


3t 


ri 


3A<pp 

Bt 


SB 

2  O 

ff  r*  at 


(45) 


(46) 


If  the  signal  is  now  integrated  with  a  passive  network,  the  detected 
has  the  form 


voltage 

(47a) 


In  fact,  the  last  term  of  this  equation  can  be  eliminated  quite  adequately 
by  using  a  second  loop  far  down  the  field,  and  subtracting  its  signal  from 
Eq.  47a,  i.  e. , 
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V 


1 

RC 


(47b) 


Now  it  is  necessary  to  relate  the  probe  signal  (Eq.  47b)  to  the  theory 
of  Section  II.  This  will  be  done  by  using  A  <p  . 

For  a  straight  magnetic  field  and  hydromagnetic  equilibrium  (see 
Ref.  2) 


p+  - -  =  constant.  (48) 

2  Ho 


In  the  present  case, 

P 


this  becomes 
(B  -  B  ) 

+  —2 - e_ 


2  Hr 


(B  +  B  )2 
o  e 


(49) 


at  the  plasma-field  boundary.  Assuming  Be  «  BQ ,  Eq.  49  shows  that 


B 

_J2. 

Bo 


/ 

1  -vl  -  (3 


1  + 


P 


(50) 


and  from  Eq.  43 


A<Pn  =  I 


Bp  dS 


=  B0  J 


1  + 


•v/i-0 


dS  , 


(51) 


where  Sp  =  the  plasma  cross-sectional  area. 

But  p  =  ,  the  transverse  energy  per  unit  volume,  and  thus 

Eq.  51  yields 


= 


2^  W± 


P  B 


o 


j  t 

Q  1 


V  dS 


2^o  W*± 


S  1  +  ^ -P 
P 


bd  i+yrj 


(52) 
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or,  by  using  Eq.  40, 


A<Pd  =  Bo  J_  (1  -/TJi)  dS  =  1Tr=  B0  (  1  -frp) 


*  rp  Bo 


1 


J3_ 

B 


o 


IT r2  B.  (1  +  A  -  A2  +  .)  . 

P  0 


(53) 


The  following  .expression  then  links  the  voltage  (Eq.  47b)  to  the  theory: 

(1  +  A  -  </a2  + '.) 


v 


IT  r2  B 

_ P  2 


RC 


®oro 


„2  2 
C3  -  *T 


r2  -  r2 
r3  rp  J 


RC  (-v/A^+I-A) 


2  2 
^3  -  *T 


r3  " 


r2 
x  o 


(^A2  +  l  -A) 


( 1+A  -  Va2^  1) .  (54) 


This  expression  gives  the  time  behavior  of  the  diamagnetic  signal  induced 
by  the  plasma  at  any  location  in  space. 
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IV.  DETERMINATION  OF  PARAMETERS  FROM  MAGNETIC  SIGNALS 

Since  six  parameters  (X  ,  r  ,  Z  ,  (Tp  +  T-)  ,  N-,  v  or  their  equiva- 

OOO  o  1  o 

lents  must  be  determined,  six  relatable  experimental  quantities  are  needed. 
As  in  the  French  work,  four  magnetic  loop  signals  are  used  to  obtain  five 
pieces  of  information: 

1.  Each  of  three  loop  signals  gives  the  maximum  value  of  the  voltage, 
vx  ,  va,  v4  . 

2.  The  difference  in  the  time  when  the  maxima  of  loop  1  and  loop  4  sig¬ 
nals  occur,  i.  e.  ,  t4  -  ti  ,  is  used 

3.  The  full  width  at  half-maximum  of  the  loop  3  signals,  plotted  as 
A3  vs  t,  which  is  denoted  At 

,i, 

where  denotes  maximum  of  quantity  and  t  denotes  time  at  which  maxi¬ 
mum  occurs. 

The  flux  tube  radius,  r  ,  is  determined  independently,  giving  a 
total  of  six  experimental  quantities.  It  should  be  pointed  out  that  these 
are  in  no  way  unique,  nor  is  there  any  reason  to  believe  at  present  that 
they  are  the  best  choice  (in  terms  of  being  most  sensitive  in  fitting  the 
data). 

Determination  of  rQ  --  There  are  several  ways  of  determining  rQ, 
such  as  electric  probe  measurements  far  down  the  flight  path,  time-inte¬ 
grated  photography,  and  the  ratio  of  magnetic  signals  from  loops  with 
different  radii  at  the  same  axial  position.  These  are  well  outlined  in  the 
reports  listed  above  (see  Ref.  1). 

The  last  method  has  been  used  in  the  present  work.  An  axial  probe 
of  radius  r2  was  inserted  into  the  vacuum  system  so  that  the  probe  and 
one  of  the  magnetic  loops  (usually  number  3  or  4)  were  coplanar.  This  is 
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shown  schematically  in  Fig.  2.  A  comparison  of  the  signals  induced  by  the 

plasma  in  these  probes  makes  it  possible  to  deduce  the  flux  tube  radius, 

r  ,  as  well  as  the  plasma  radius  r  .  This  method  is  described  below: 
o  P 

Equation  47b  gives  the  voltage  recorded  by  a  single-turn  magnetic 
loop  surrounding  the  vacuum  system.  When  the  axial  probe  is  used,  two 
cases  are  to  be  considered:  ry  >  rp  and  ry  <  rp .  Clearly,  case  1  is 
covered  by  Eq.  47b  where  rx  is  replaced  by  r?  .  It  is  equally  clear  that 
the  ratio  of  the  recorded  voltages  yields  no  information  about  rp  (and  thus 
about  rQ). 

.For  case  2,  however,  the  flux  through  the  axial  probe  as  a  function 
of  time  is  given  by 


<p(t)  =  TT  rl  (Bq  -  Bp)  , 
from  which  it  follows  that 


v 


axial 


RC 


TT  r|  B 
_ o 

RC 


(55) 


This  expression  is  dependent  on  rp  . 
voltage,  vi00p’  the  Pr°be  voltage, 
it  is  found  that 


From  the  ratio,  p  ,  of  the  loop 

v  .  measured  at  the  same  time, 
axial 


p  r  l  r  l 

(r32  -  r?  +  p  r| ) 


(56) 


and  that 


v  RC 

-Lqpjb _ 


trB 


o 


(57) 
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The  question  arises  whether  the  axial  probe  disturbs  the  plasma  mark¬ 
edly  in  this  measurement.  This  was  investigated  both  in  the  French  wo~k 
and  in  the  present  work.  In  the  former  investigation,  .r  as  measured  by 
the  above  techniques  was  approximately  10fo  smaller  than  that  found  by 
other  methods.  In  the  present  case,  the  magnetic  loop  signals  were  es¬ 
sentially  unchanged  by  the  presence  of  the  probe;  nevertheless,  the  four 
loop  signals  used  in  this  analysis  were  recorded  with  this  probe  withdrawn 
and  a  separate  run  was  made  for  the  purpose  of  measuring  r  .  (see  Fig. 

3). 

,  Extraction  of  Data  from  Magnetic  Loop  Signals- -  Five  quantities  re¬ 
main  to  be  obtained  from  the  magnetic  loop  signal:  the  maximum  voltages 
from  loops  1,  2,  and  4  (vx  ,  v2  ,  v4  ),  t4  -  U  ,  and  the  full  width  at  half¬ 
maximum  of  loop  3  signal  plotted  vs  t. 

The  first  four  of  these  items  may  be  read  directly  from  oscilloscope 
tracings.  The  last  quantity  is  slightly  more  cumbersome  to  obtain;  the 
oscilloscope  trace  (v3  vs  t)  must  be  converted  into  a  plot  of  A3  vs  t. 

Using  Eqs.  40  and  54,  it  can  be  shown  that 


where 


(58) 


(59) 
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These  equations  permit  the  data  (in  the  form  of  v3  vs  t)  to  be  plotted  as 
A3  vs  t  from  which  A  t  [-^J  can  be  determined  graphically.  Ax ,  A2  , 
and  A4‘  (which  correspond  to  vx' ,  v2 '  and  v4’  )  can  also  be  calculated. 

The  correspondence  is  due  to  the  fact  that  v  is  a  monotonic  increasing 
function  of  A. 

The  information  necessary  to  carry  out  the  analysis  has  now  been 
obtained.  It  should  be  mentioned  that  the  calculations  described  in  this 
section  (Eqs.  56,  57,  58,  and  59)  were  done  on  a  computer  as  was  a 
large  fraction  of  the  work  to  be  discussed  in  the  following  section.  These 
computer  programs  are  described  in  the  appendices. 
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V.  ANALYSIS 


The  analysis  is  carried  out  essentially  by  graphical  methods,  in 
which  successive  approximations  are  used  to  give  the  best  fit  to  the  data. 

The  analysis  procedure  will  be  outlined  first  and  is  based  on  the 
assumption  that  the  necessary  graphs  are  available.  Then  the  means  of 
obtaining  these  curves  will  be  discussed.  The  required  graphs  are  as 
follows: 

1.  The  family  of  curves  f"  {a,  Z1)  =  constant  wnere  f'  {a,  z')  is  the 
maximum  of  the  function  (Eq.  35)  with  respect  to  T. 

2.  The  family  of  curves  T  r  =  constant  as  a  function  of  {a,  Z1)  where 
T  ’  is  the  reduced  time  at  which  the  maximum  of  f  (  A>  a,  Z1)  occurs,  in¬ 
dependent  of  the  value  of  f  for  the  particular  {a,  Z1,  A  )  . 

3.  The  family  of  curves  At^  =  constant  as  a  function  of  a,  z'  ,  where 
AT|  is  the  full  width  at  half-maximum  of  the  function  f  (a,  z',  A),  plotted 
as  a  function  of  T  . 


The  following  ratios  are  constructed  next: 


Ql2 


(60a) 


and 


Ql3  - 

Ks  = 

z4  -  zx 

f*  {a,  Za'  ) 

Z4  -  Z 

^■2 

f*  (a-,  Z{) 

Ar 

II 

f*  (a,  z[) 

A4:< 

f:;:  ia-,  Zx') 

f— 1 
< 

1 

(60b) 


(60c) 


(60d) 


The  equality  of  the  ratios  can  be  seen  from  Eqs.  30  and  37. 
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It  is  noteworthy  that  the  first  terms  to  the  right  of  the  equal  signs  are 
constructed  of  variables  from  the  theory  while  the  second  terms  are  made 
up  of  experimental  determinable  quantities.  (The  values  are  calculated 

from  the  values  of  v-  ;  Z-x  is  the  distance  from  the  center  of  the  gun  coil 
to  the  loop  positions).  These  are  experimental  values  for  the  ratios  in 

i 

Eq.  60  for  which  a  best  fit  in  the  (q-,  Z  )  plane  must  be  found. 

A  starting  set  (q-,  Zx )  is  chosen  arbitrarily.  The  corresponding 
f'1'  (q,  Zx  )  is  found  from  the  graphs  of  f "  =  constant.  Then  a  value  of 


f"  (O  ,  Z*  )  =  K4GXP  f*  (Q;  Zj  ) 


is  calculated.  (Superscript  "exp"  indicates  experimentally  determined 

*  1  1 

values.  )  With  this  value  of  f  (co,  Z4  )  the  coordinate  Z4  which  corresponds 
to  the  value  [<?•,  fv  (o,  Z4  )  ]  is  located  on  the  graph.  Since  from  Eq.  60a 
Z2'  =  z[  +  Q12  (Z 4  -  z{  ),  it  is  possible  to  find  f (a>  Z2  )  at  the  coordi¬ 
nates  (as  Z2  )  and  construct 


r‘  (a-,  Z2'  ) 

^  =  f;:"  (o,  Zj*  ) 

based  on  an  estimate  of  (a ,  Zl  )  and  compare  it  to  K2  ‘  r  derived  from  the 
data.  This  procedure  is  repeated  and  the  choice  of  (cr,  z|  )  is  adjusted 

0XD 

until  good  agreement  is  achieved  between  K2  and  K2  *  .In  this  manner 
the  range  of  ( a z[  )  is  limited. 

These  values  of  (a-,  Z-'  )  are  used  with  the  second  set  of  graphs,  T  = 
constant  as  a  function  of  (a  ,  Z1  ).  The  experimental  information  to  be  used 
is  the  difference  in  the  times  at  which  the  maxima  v4  and  vx’  occur,  i.  e.  , 
t4  -  ti'  .  From  Eq.  21b 


(61) 
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>|t  |  I 

Locating  t4  at  (a,  Z4  )  and  Tj  at  (a,  Zl  ),  £  may  be  evaluated. 

It  is  now  possible  to  limit  further  the  allowable  values  of  (a ,  Z?)  by 
using  the  full  width  at  half-maximum  of  the  signal  from  loop  3.  A  at 
(a ,  zl  )  may  be  found  from  the  curves  A  Tx  =  constant,  where  Eq.  60b 
becomes 

Z3  =  Z\  ±  Q13  (Z4  -  Z\  ) 

1 

For  the  proper  values  of  {a,  Z±  ) 


>;c 


r  i 

1  1  1 

"_A_‘ 

t> 

s 

<1 

I 

II 

L  *J 

3  i 

2 

Generally  it  is  necessary  to  repeat  these  three  steps  several  times 
in  order  to  achieve  the  best  fit  to  the  data.  This  sequential  processing  per¬ 
mits  one  to  limit  to  a  fairly  narrow  domain  the  allowable  range  into  which 
the  characterizing  parameters  fall.  ( Waelbroeck  et  al,  Ref.  1,  find  ac¬ 
ceptable  fits  with  X  ±  8%,  Ni  ±  6%  ,  Z0  ±1.5  cm,  (T  +  Tj  ±  15%  , 
v0  ±  5%). 

As  a  final  step  these  parameters  are  used  to  calculate  the  voltages 
from  Eq.  54  and  compare  them  to  the  traces  recorded  by  the  oscilloscope  . 


I 
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VI.  COMPUTER  ANALYSIS 

The  computer  analysis  to  be  described  below  was  carried  out  on  the 
IBM  7094  computer  at  APL/JHU.  It  is  by  no  means  the  most  sophisticated 
means  of  conducting  the  analysis  but  is  designed  to  eliminate  some  of  the 
mathematical  drudgery  which  the  method  described  in  Section  V  requires. 

The  programming  is  divided  into  five  parts: 

1.  The  calculations  of  r  ,  rQ,  A3  vs  t,  ,  A2 A4'  and  the  ratios 
K2 1  K4  #  Q12  i  Q13 »  (Subsec.  V). 

2.  The  calculation  of  the  family  of  curves  fv  (a,  z')  =  constant. 

3.  The  calculation  of  the  family  of  curves  T  *  (os  z')  =  constant. 

4.  The  calculation  of  the  family  of  curves  A (0,  Z1)  =  constant. 

5.  The  calculation  of  the  voltages  using  the  parameters  from  the  anal¬ 

ysis,  the  local  velocity  and  the  volume  density  m. 

Each  of  these  calculations  will  be  treated  in  turn. 

The  Fortran  listings  of  these  programs  are  given  in  Appendix  V. 


-  29  - 


Th»  Johni  Hopklnt  Univ«riliy 

applied  PMvaica  LAaoRAToay 

Sllv*f  Spring,  Maryland 


Reduction  of  Magnetic  Loop  Data-  -  This  program  is  the  most  straight¬ 
forward  of  the  five.  The  mathematics  is  outlined  in  Subsec.  .IV  and  needs 
little  further  explanation.  In  this  program,  the  Mks  system  of  units  is  used 
and  the  indicated  times  are  in  microseconds. 


The  input  data  necessary  to  this  program  are: 

1.  ;  VOUT  (the  peak  voltage  from  the  magnetic  loop  for  measure- 
of  r0) 

2.  V  ;  VIN  (the  peak  voltage  from  the  axial  probe  for  measurements 

axiai 

of  ro) 

3.  r2  ;  Rn  (radius  of  the  axial  probe;  Vn  ^  is  used  for  more  than  one 
turn) 

4.  R,  R* ;  R,  RPR  (the  resistance  in  integrator  from  the  loop  and  the 
axial  probe,  respectively) 

5.  C,  C* ;  C,  CPR  (the  capacitance  in  integrator  from  the  loop  and  the 
axial  probe,  respectively) 

6.  r3;  R3  (flux  conserver  radius) 

7.  rx;  R1  (magnetic  loop  radius) 

8.  V]/'5,  v2'\  v3v,  v4";  VMAX  1  (voltage  maxima  from  loops) 

9.  t2  ,  t3  ,  t4‘;  TMAX  1  (time  at  which  the  voltage  maxima  occur) 

10.  Zlt  Z2,  Z3,  Z4  ;  Z1  (location  of  loops  from  the  center  of  the  gun 
coil) 

11.  B0  ;  BO  (magnetic  field  intensity  of  drift  field) 

12.  v3  vs  t  (loop  3  voltage  and  the  corresponding  times) 

The  output  data  are  shown  in  the  appendix  for  a  typical  case.  They  are: 

(1)  r0 ,  (2)  r  ,  (3)  A1‘P,  As",  A3",  A4'',  (4)  A3  vs  t,  (5)  K2,  K4,  Q12 ,  Q13 
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Calculation  of  Curves  f  '~  (a,  Z  )  =  Constant--  The  equal- value  curves 
of  f'  (a ,  Z 1  )  are  obtained  in  the  following  way: 


i  -  B  /3  7*  -  rvT  ~ 

f  (  A,  a,  Z  )  =  A  exp  -  - — -  (Eq.  35). 

df 

In  order  to  find  the  maximum  of  f  with  respect  to  T ,  —  is  set  equal  to 
zero  yielding 

«  z  -  aT  dA  Z1  -  aT  5  dA  . 

2  A  dr  +  2a  A  "  7  ~&T  ‘  0  '  <63) 


Solving  Eq.  63  as  a  quadratic 


•  •  Z1  -  or 

10  in  L  A  J  ' 


and  requiring  that  this  argument 


remain  finite  as  a  -  «  ,  i.  e.  ,  [(T  +  T.)0  -  o]  gives 

e  1 


z'  -  aT* 


-  a  +  Va  +  “T  ^1  -  A 


r  *  2/3U 

2  [i  -  a  r 


(at  maximum)  , 


where 


1  -A 


2/3]  £ 


(Eq.  24) 


Eq.  35  may  be  written  in  the  form 


Z  <y  T  „  ,  5  ,  A 

- - -  =  In  f  +  —  In  A 

A  J  3 


Equations  64  and  65,  at  the  maximum,  are  combined  to  give 

5  ...  5  *  ...  -2/3 

7  +  In  r  +  7  In  A  1  -  A" 

b  0  J  . 


-  In  f"  -  |  In  A"' 


The  curves  f*  (a,  z'  )  =  constant  may  be  found  using  this  expression. 
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First  the  constant,  which  may  range  between  0  and  1  is  chosen;  then 
a  set  of  values  of  a  are  generated  (using  Eq.  66)  for  a  series  of  values 
A*  greater  than  1.  The  corresponding  Z1  is  then  calculated  from  Eq.  64 
and  the  expression  t(  A)  is  given  by 


...  2/3 

,  2/3 

T  "  = 

A"  +  2 

A”'  -  1 

(Eq.  23)  . 


;'c  :!c 

The  restrictions  on  the  allowable  values  for  the  pairs  (f ",  A  )  are 
readily  found  from  Eq.  66  to  be 


-  e/s  ,  s/3 

e  s  f*  A"  <  1  .  (67) 


This  domain  is  shown  graphically  in  Fig.  4. 

Since  the  equal- value  curves,  f'  (a,  Z  )  =  constant  and  t"  (q,  z')  = 
constant,  were  computed  in  the  same  program,  discussion  of  this  program's 
input-output  data  will  be  deferred  until  the  end  of  the  next  subsection. 


W  | 

Calculation  of  the  Curves  T  '  (a,  Z  )  =  Constant--  This  calculation 
is  conducted  in  a  manner  similar  to  that  of  the  preceding  section.  First, 
a  value  A  =  constant  is  chosen  and  this  implies  a  value  for  T  =  constant 
(Eq.  23);  a  series  of  values  of  a  are  then  introduced  into  Eq.  64  to  find 
Z1.  Thus  the  sets  (»-,  Z1)  for  which  T ‘‘  =  constant  are  generated  and  this 
process  is  repeated  for  successive  values  of  A". 

There  are  two  parts  to  the  input  data  to  this  program:  input  data  for 
f  ’  (a,  Z1)  =  constant  and  input  data  for  T  "  (o-,Z'  )  =  constant. 
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A.  The  input  data  for  f'  (a-,  z')  are: 

1.  Initial  value  of  f  :  FMAX 

2.  Initial  value  of  A*  :  XLAM 

3.  Incremental  change  in  fv,  (Af*)  :  DELFMX 

4.  Incremental  change  in  A*  (AA  )  :  DELAM 

5.  M,  the  number  of  values  of  f’,:  to  be  used. 

6.  N,  the  number  of  values  of  A  to  be  used. 

B.  The  input  data  for  T*  (a*  Z1  )  =  constant  are: 

1.  Initial  value  of  a  :  XALPA 

2.  Initial  values  of  A  *:  YLAM 

3.  Incremental  change  in  a, (Ace)  :  DELALP 

4.  Incremental  change  in  A?  ( AA*)  :  DEYLAM 

5.  L ,  the  number  of  values  of  A  to  be  used. 

6.  K,  the  number  of  values  of  a  to  be  used. 

Similarly,  the  output  data  are  divided  into  two  parts. 

A.  The  output  data  for  f‘  (a,  z)  are: 

1.  The  value  of  f" 

|  jJ(  jjj  jij 

2.  A  listing  of  ce,  Z  ,  t  ,  A  corresponding  to  this  value  of  f 

B.  The  output  data  for  (g»  z')  =  constant  are: 

1.  The  value  of  T* 

2.  A  listing  of  a  and  Z1  corresponding  to  this  value  of  T * 

For  the  purposes  of  the  analysis  these  data  were  graphed.  Typical 
curves  are  shown  in  Figs.  3  and  5. 


-  34  - 


V  -f?.. 


■nprra 


2  (33.0  CM) 


TIME  ( 


TK*  John*  Hopkini  Univ«filty 

APPLIED  PHYSICS  LABORATORY 

Silver  Spring.  Maryland 


Calculation  of  the  Curves  AL  =  Constant--  This  is  by  far  the  most 
complicated  calculations  and  probably  the  most  interesting  from  the  pro¬ 
grammer's  point  of  view. 

The  full  width  at  half-maximum,  AT^,  is  given  by 


Ar^  =  ra  -  7i  ,  (68) 

where  r2  and  Tx  are  the  reduced  times  at  which  the  function  f  (  A  ,  a-,  z') 
achieves  half-maximum  values,  £  f The  widths  are  found  in  the  following 

way:  A  value  of  A"  (and  thus  t',c),  and  a  value  of  a  are  chosen;  Z1  is  cal- 

*  i  * 

culated  from  Eq.  64.  Next  f  (A”",  q,Z  )  =  f  is  computed.  Finally  the 
roots  of  the  following  equations  are  found: 


(69a) 


with 


and 


1  <  Ax  £  A  ' 

-5 /  3 

=  As  exp  - 


(69b) 


with 


A"  £  Ag  s  M  A  ’  , 


where  M  is  some  sufficiently  large  multiplier.  Having  found  Al  and  A2, 
T i,  rs,  and  A  Ti  are  found  directly.  In  the  present  case,  a  bisectional 
method  was  used  to  find  the  roots  Ax  and  As  of  Eqs.  69a  and  69b  to  an 
accuracy  of  0.  001.  (For  this  purpose,  t1  and  t2  were  expressed  in  terms 
of  Al  and  Az  according  to  Eq.  24. 
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For  each  value  of  A  ,  a  list  of  a,  Z*  and  A  T:  is  calculated,  the 
values  of  a  used  in  this  program  were  the  same  for  each  A  (  A  ranged 
from  1  to  10). 

i 

The  graphs  A  (o  ,  Z  )  =  constant  are  obtained  by  first  plotting  the 
family  of  curves  a  =  constant  as  a  function  of  AT:  and  Z*,  and  then 
interpolating  to  find  the  curves  desired.  These  curves  are  presented  in 
Fig.  6. 

The  input  data  for  this  program  consist  of: 

1.  The  accuracy  to  which  the  roots  A ^  and  A2  are  to  be  found:  EPS 

2.  The  initial  value  of  a  :  XALFA 

3.  Initial  values  of  A  '  :  XLAM 

4.  Incremental  change  in  a,{ Aa)  :  DEALF 

5.  Incremental  change  in  A  ,(A  A  ')  :  DEXLM 

6.  M,  the  number  of  values  of  A  to  be  used. 

7.  N,  the  number  of  values  of  a  to  be  used. 

The  output  data  for  this  program  consist  of: 

1.  A*  :  LAMAX 

2.  o  :  ALPHA 

3.  Z'  :  ZPRME 

4.  ATi  :  WIDTH 

T? 

5.  Calculated  value  of  ^f'  at  root  A;  :  FI 

6.  Calculated  value  of  |f'  at  root  A2  :  F2 

7.  Ai  near  root  found  by  the  program:  XLAM1 

8.  A2  for  root  found  by  the  program:  XLAM2 

9.  r  :  FMAX 
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Reconstruction  of  Magnetic  Signals- - When  the  curves  described  above 
are  available,  the  analysis  of  Section  V  can  be  carried  out  to  find  the  five 
remaining  unknown  parameters  which  describe  the  plasma.  The  final  step 
is  to  reconstruct  the  four  magnetic-loop  signals  based  on  the  theory  and 
using  the  extracted  parameters. 

The  parameters  which  are  found  from  the  analysis  are  \0  ,  a, 
r0  ,  N.,  and  Z0  .  The  loop  signals  are  reconstructed  by  substituting  the 
above  data  in  Eq.  37  to  find  A  as  a  function  of  t,  by  first  calculating  A(A) 
and  then  converting  A  to  the  time  domain.  Using  these  results,  v  as  a 
function  of  time  is  then  calculated  from  Eq.  54. 

In  the  program  used  here,  the  volume  density,  n.(Z,  t)  and  the  local 
velocity,  v(Z,  t),  were  also  calculated. 

The  input  data  for  this  program  consist  of: 

1.  Value  of  a  :  ALPHA 

2.  Initial  value  of  A  :  XLAM 

3.  Incremental  steps  in  A  :  DX 

4.  Value  of  X0  :  XLAMO 

5.  Value  of  Z0  :  ZO 

6.  N,  the  number  of  steps  in  A  to  be  made 

7 .  Val  ue  of  r0  :  RO 

8.  Val  ue  of  rx  :  R1 

9 .  Val  ue  or  r3  :  R3 

10.  Value  of  Bq  :  BO 

11.  Integrator  time  constant,  RC:  RC 

12.  Cent  er-of-mass  velocity,  vx  :  VO 

13.  Ion  mass,  m.:  XM 

l 
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14.  Total  number  of  ions,  N^:  XN 

15.  Value  of  £  :  XI 

16.  Loop  location,  Z,  and  the  last  loop  location:  Z,  ZLAST. 

The  output  data  are: 

1.  The  characteristic  parameters 

2.  The  voltage  V,  the  time  t.  A,  T,  A,  density  n^ ,  and  velocity  v. 

A  separate  listing  is  given  for  each  loop  location,  Z. 

This  now  makes  the  calculated  signals  available  to  be  compared  to 
the  oscilloscope  traces. 
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VII.  RESULTS 


The  analyses  of  two  plasma  puffs  are  presented  here.  The  initial 
conditions  in  the  two  cases  were  significantly  different,  these  differences 
being  reflected  in  the  final  plasma  characteristics. 

Figure  4  shows  a  typical  set  of  magnetic  loop  signals  to  which  the 
foregoing  analysis  was  applied.  The  positions  of  the  loops  relative  to  the 
gun  face  are  indicated. 


Figure  5  shows  the  fit  obtained  for  one  case.  The  solid  curves  rep¬ 
resent  the  loop  signals  while  the  points  indicate  the  values  computed  by 
the  theoretical  method.  This  plasma  can  be  characterized  as  dense  and 
cool.  As  can  be  seen  from  the  figure,  the  fit  improves  as  the  plasma  ad¬ 
vances  through  the  drift  field.  The  plasma  characteristics  arrived  at 
are: 


and 


(TP  +  T:)  =  6  x  103  °K 

1  o 

v0  =  3.45xl06  cm/sec 
IT  =  2.  7  x  1019  ions 
X0  =  2,16  cm 

Z0  =  17.1  cm 


Of  some  concern  is  the  fact  that  the  computed  curve  decays  more 
rapidly  than  the  signal  itself.  No  satisfactory  explanation  for  this  has  been 
found.  A  possible  explanation  may  be  that  the  field  from  the  gun  influences 
the  magnitude  of  the  first  loop  signal  significantly,  in  which  case  the  value 
of  A*  is  different  from  the  value  whicn  wuuid  be  obtained  otherwise.  The 
resulting  "best"  set  of  parameters  to  fit  the  data,  especially  for  a  low- 
temperature  plasma,  are  adversely  affected. 
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The  low  temperature  and  high  density  calculated  for  this  case  sug¬ 
gest  the  possibility  that  the  assumptions  for  the  computation  may  not  be 
fulfilled.  Initial  trapped  field,  recombination,  or  diffusion  may  be  play¬ 
ing  a  role.  Further  experimental  measurements  are  expected  to  clarify 
this  picture. 


Figure  6  shows  the  fit  to  the  second  plasma  puff.  The  fit  to  these 
data  clearly  is  considerably  better.  This  plasma  has  the  following  char¬ 
acteristics: 


(T  +  T.) 
e  1  0 

v0 


2 . 2  x  1 04  °K 
7.  5  x  106  cm/sec 


N.  =  1 . 07  x  ID19  ions 

i 

X0  =  4.45  cm 

Z0  =  0.  64  cm 


The  fit  to  all  four  signals  is  quite  good,  with  a  marked  improvement 
in  the  fit  to  the  trailing  edge  of  the  pulses.  In  support  of  the  previous  dis¬ 
cussion,  it  can  be  seen  that  better  overall  fits  were  obtained  by  sacrificing 
the  fit  to  the  magnitude  of  the  first  signal. 


It  should  be  noted  that  the  0  calculated  for  this  case  comes  out  to 
0.  86,  which  though  physically  reasonable  and  not  in  disagreement  with 
the  observed  puff  size,  does  cast  some  doubt  on  the  numbers  above. 
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VIII.  DISCUSSION 

Further  analysis  of  the  plasmas  derived  from  the  gun  is  planned, 
both  with  the  longitudinal  field  and  the  magnetic  octupole  field  using  this 
model. 

The  shifting  of  the  pick-up  loops  downstream  will  be  attempted  to 
estr  blish  whether  the  gun  influences  the  signal  from  the  first  loop  suffi¬ 
ciently  to  affect  the  final  results. 

It  is  noteworthy  that  although  the  gun  and  the  drift  field  used  here 
were  quite  similar  to  the  arrangement  used  by  the  French  investigators, 
the  plasmas  produced  here  are  typically  denser  and  of  lower  temperature 
than  those  produced  by  the  prototype  gun.  It  is  quite  likely  that  these  dis¬ 
crepancies  are  due  to  the  difference  in  gun  frequency  (the  present  gun 
has  a  ring  frequency  of  80  kc  while  the  French  gun's  frequency  was  a 
factor  of  two  higher. 


-  43  - 


The  Johru  Hopkini  Univer*ity 
APPLIED  .PMYtICt  LADORATORY 

Silver  Spring.  Maryland 


APPENDIX  I 


Solution  of  Equation  1  Ob 


The  following  relationship  was  found: 


n.V  ■  -  X*  J 


2Z 

-  n-  "  T 

3  1  X 


Hi-1 


LX 


dZ  , 


(10b) 


or 


X  N. 

n.  V  =  -  - -  r 

1  fW  J 


\3  / 

\  j 

2Z2  / 

Z  \  1  / 

z  \ 

exp  -l 

r/  ■  exp 

X4  ' 

x  /  x2  \ 

A  / 

dZ  . 


Let 


Z2  =  x,  then  dz  = 


dx 

2^T 


Thus, 


n.  V  = 
i 


X  N 

_ j 

Vir 


J 


A 

X4 


exp  - 


(?) 


2  X‘ 


exp  - 


(?) 


dx 


(10c) 


Integrating  the  first  term  on  the  right  yields 


J 


X 

H  ^ 

X 

1  „  l 

r  “ 

V 

exp' 

_f 

dx  =  -  ^2  exP  " 

X2 

+  2X5  J  77 exp ' 

A 

F 

dx  +  c  , 


where  C  is  the  constant  of  integration. 
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Equation  10c  now  becomes 

X'n 

n.  V  =  - —  Z  exp 

1  x'yi 

Since  V  is  the  velocity  relative  to  the  center  of  mass  of  the  plasma 
puff,  the  velocity  of  the  particles  at  a  distance  +z  must  be  equal  to  the 
velocity  of  particles  at  -z.  Thus,  the  constant  C  must  be  identically 
zero. 

i 

7  X 

V  =  — —  Q.  E.  D.  (11) 


_Z 

X 


+  c 


(10d) 
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APPENDIX  II 


Solution  of  Equation  12 


„  n .  m .  Z  X 

|£.  +  -L-i -  =  0 

SZ  X 


r  n.m .  ZX 

J-i-I - 

X 


dZ  = 


Na  X 
A 


f  - 


It  should  be  noted  that  taking  N  outside  the  integral  sign  is  again  invoking 
the  assumption  that  r^  is  a  slowly  varying  function  z. 

Again  letting  Z2  =  x, 


V 


n  1  X  , 

m.  f  — V  exP  "  dx  ' 

i  J  2XS  X 


m.NA 
l  A  %  " 

-  a  exp 

2^r 


Z_ 

"  .  X. 


+  c  , 


Z  -  ±  03  ,  p  -*  O  . 


Tnus, 


m.NA 
l  A 

2 


X  exp  - 


_Z_ 
_  X_ 


Using  Eq.  8,  we  have 


p  =  ~  n.m.  X X 

2  i  i 


Q,  E.  D. 
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APPENDIX  III 


Let 


Then 


where 


Solution  of  Equation  22 


3  A 


cl2  A 
dr2 


P 


dA 

dr 


1  . 


dp  _  d2 A  =  dp  d A 
dr  dr2  dA  dr 


d2  A 
d  t2 


P  P 


i 

P 


d£_ 

dA 


(22) 


Equation  22  now  may  be  rewritten  as  follows: 


Multiplying  both 


3p'p  A  +  p2  =  1  . 

n  -  1 

sides  of  Eq.  22a  by  A 

O  a  n  1  ,2  An_1 

3  A  p  p  +  p  ^  = 


yields 

A"-1  . 


(22a) 


(22b) 


But  (Cp2An)  =  2  Cp  p' An  +  n  Gp2  An_1  , 

where  C  is  some  constant.  This  expression  is  exactly  equal  to  the  left 
side  of  Eq.  22b  if 
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2C  =  3  and  nC  =  1 

C  =  3/2  n  =  2/3 

Equation  22b  may  be  written  as 

~t  (3/2  p2  As/3)  -  A  1/3  • 
d  A 

Integration  of  this  equation  gives 

3/2  p2  A2/3  =  3/2  A 2/3  +  L 
or 

p2  =  1  +  M  A"2/3  , 


where  M  is  the  constant  of  the  integration. 


dA 

P  =  dr  = 

(1  +  M  A’2/3)  ^ 

d  T  = 

d  A 

(1+  M  A"2/3)^ 

Let 

V 

% 

II 

> 

Then 

d  A  = 

3y3  dy  , 

and 

,  _  3  y2  d  y 

3  y3  d  y 

dT  -  1 

(1  +  My"2  r 

(y2  +  M)^ 

r  ~  (y3  - 

•  2M)  (y2  +  M)  ^  , 

T  =  (A2/ 

3  -  2M)  (A3/s+  M)^  , 

(22c) 
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at  t=r  =  0  A  =  1  . 


.  ‘  .  M  =  -1  or  1/2  . 

Using  Eq.  22c,  Eq.  25  of  the  text  becomes 


(T  +  T.) 
e  i 


(T  +  T.) 
e  l 

o 


and  since  ratio  must  be  positive,  M 


Thus  we  *<nd 


T 


(A2/3  +  2)  ( A3^3  -  1)^  , 


dA 

dT 


(  1  -  A -a/3)*  . 


Q.  E.  D. 


(23) 

(24) 
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APPENDIX  IV 

The  Compensated  Magnetic  Loops 

The  compensated  loops  used  by  the  French  investigators  effectively 
consisted  of  two  loops  of  equal  areas  in  opposition.  One  of  these  loops  en¬ 
circles  the  plasma  while  the  other  encloses  an  area  outside  the  plasma  but 
inside  the  flux  conservcr.  (see  sketch  below) 


Sch  matic  of  Compensated  Magnetic  Loop  Arrangement 


The  flux  passing  through  the  loop  may  be  calculated  by  means  of 

<p(t)  =  r 2  (B  -  B  )  +  n(  rl  -  r2)  (B  +  B)-T7r2(B  +  B) 
pop  poe  eo 

where  B  =  field  intensity  from  solenoid 
o 

B  =  field  generated  internal  to  plasma 
P 

and 

B  =  field  due  to  flux  excluded  by  plasma 
c 

The  first  term  above  represents  the  flux  through  the  plasma  cross 
section.  The  second  term  is  the  flux  in  the  region  between  the  plasma  and 
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the  loop  which  encircles  the  plasma,  while  the  last  teim  is  the  flux  encircled 
by  the  compensating  element  whose  effective  area  is  equal  n  rf  . 


<p  (t)  = 


-  nr2  B  -nr2  B 
P  P  pc 


But 


nr  B  =  A  <p  =  rr  (  m  -  r 2  )  B 
P  P  P  pc 


and  thus 


cp  (t)  =  -  A<p 


1 


2  2 

"3  -  "p  J 


The  e.  m.  f.  ,  £  ,  is  given  by 


^  9  t 


1  + 


r2  -  r  2 
r3  rp 


5Aipf 

St 


and  the  integrated  signal 


RC 


1  + 


2 

r3 


P 


A  <p 


This  type  of  magnetic  loop  has  the  obvious  advantage  of  eliminating  the 

effect  of  the  B  field. 

o 
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APPENDIX  V 


Computer  Listings 


52 


REDUCTION  or  SINGLE  LOOP  DATA’RO  CALCULATION 
01 n  MS  ION  VII56)  *  r  3 ( 56 ) , A  3  X ( 5  6 ) *  03X156) 

DXXMS,T,U,X,Y)  =  (((S/n**?.)-((U/n**2.)MX/Y))/(  (  l.MX/Y  )  )  *(  (S 

1/ T ) **2. I-( ( U/T ) **2. )  ) 

AXXM  X  )  =  (  1  {  X**2.  )  )  /  (  2.  *X  ) 

55  PfcAO  INPUT  TAPE  5 , 30  0  ,  VOU  T  ,  V  I  N  ,  R  2  ,  R  P  R  ,  CPR 
READ  INPUT  T APE5 ,300  ,R3,R1 ,DO,R ,C 
RtAU  INPUT  TAPE  5 , 30 5 , VMA X  I , V  MAX  2 , VM AX 3 , VMAX4 
READ  INPUT  TAPE  5 , 30  5  ,  TMAX  1  ,  T  MAX  ?  ,  T  MA X 3  ,  TM AX4 
READ  INPUT  TAPE  5, 305 ,Z  I , 12 , 23 ,74 
PI  =  3. 1415927 
RATI 0= VOUT /V  IN 

RPSO=(R3**2.  )*(R2**2.)*(RATIO*R*C)/(  ( RPR  *C  PRIM  (R3**2.  )-(Rl**2.)+( 
IR2**2. ) *(  (RAT I  0* R*C  )  /(PPR*CPR ) )  )  ) 

RO  =  SQRTF  (  RPSQ- (  (  (R*C*VOUD  M  (  R  3**2 .  )  -  (  RPSQ  II)  /  (  P  l  *B0  ♦  (  (  R  3**2 .  )  - 1  R1 
1**2. ) )  )  )  ) 

RP=SORTF(RPSO) 

RETA=1.0-( (RO/RP )**4. ) 

WRITE  OUTPUT  TAPE  6 , 308 , RO , R P , BO , BE T A 
XK  =  ( P I  * ( R0**2  .  )  *R0  )  / ( R  *C ) 

DMX 1 =  DXXF(R3,R0,R1,VMAX1,XK) 

DMX2=  OXXF ( R 3, RO,R l  ,  VMAX2, XK ) 

OM  X  3  =  DXXF ( R  3 , R  0  ,  R  1  , VMAX3  » XK  I 
D  M  X  4  =  DXXF(R3,R0,R1  ,VMAX4,XK) 

AMAX 1 =  AXXF ( DMX  1  ) 

AM AX2  =  AXXFIDMX2J 
AM AX  3=  AX  X F ( DMX  3  ) 

AMAX4=  AXXF(DMX4) 

WRITE  OUTPUT  TAPE  6 , 30 1 , AM AX  1 , A M AX 2 , AM AX3 , AM AX4 , VHAX 1 , VMAX 2 , VMA X 3 , 
1 VMAX4 , TMAX1 , TMAX 2, TMAX 3 , TMAX 4 
READ  INPUT  TAPE  5,302,N 
READ  INPUT  TAPE  5 , 30 3 , ( V 3 ( I ) , I = 1 , N ) 

READ  INPUT  TAPE  5 , 3 0 3 , ( T 3 ( I ) , I = 1 , N ) 

WRITE  OUTPUT  TAPE  6,304 
DO  5  I = 1 , N 

D3X( I )=DXXF(R3,R0,R1 ,V3( I ) ,XK) 

5  A  3  X (  I  )=AXXF(D3X(  I  >  ) 

WRITE  OUTPUT  TAPE  6,306,(1  ,A3X(I  ),T3(I  l,V3(  I 
Q 1 2  =(22-Zl)/(Z4-Zl) 

013  = ( Z  3- Z 1 ) / ( Z4-Z  1 ) 

YK  2  =  ( AMAX2 ) / ( AMAX 1 ) 

YK4  =  ( AMAX4 ) / ( AMAX 1 ) 

WRITE  OUTPUT  TAPE  6 , 307 , Q1 2 , 0 1 3 , YK2 , YK4 
GO  TO  55 

300  FORMAT  (5F10.6) 

301  FORMAT  ( 7X , 6HAMAX 1 =F7 . 3 , 3 X , 6HAMA X2  =  F7 . 3 , 3X , 6H A MAX 3=F 7. 3 , 3X , 6  HA 

1MAX4=F7. 3/ ( 7X,6HVMAX1  =  F7.3 , 3 X , 6H VMAX 2= F7 . 3 , 3X , 6HVM AX  3= F7 . 3 , 3 X , 6HVM 
2AX4=F7. 3) / {7X,6HTMAX1=F7.3,3X,6HTMAX2=F7.3,3X*  6HTHAX  3=F7.3,3X,6HTM 
3AX4=F7.3) ) 

302  FORMAT  ( 5X,  12  ) 

303  FORMAT ( 7F 1 0 . 6 ) 

304  FORMAT  ( 2 X , 1 H I  ,7X,2HA3»8X,2HT3,8X,2HV3,//) 

305  FORMAT  (4F10.6) 

306  FORMAT  (  IX , I  2 , 3X , F 8. 3 , 2X , F 8 . 3 , 2 X , F8 . 3 ) 

307  FORMAT  ( 7X , 4H0 1 2=F 7 . 3 , 3X , 4HQ 1 3  =  F 7 . 3/ ( 7X, 3HK2  =  F 7. 3 v 3X  , 3HK4=F 7 . 3 1 ) 

308  FORMAT { l HI , 6X , 3HR0=F 10. 6, 6HME TER S , 3X , 3HRP= F 1 0. 6 , 6HME TERS , 3X , 3HB0=F 
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RE DUCT  ION  OF  SINGLE  LOOP  DATA'RO  CALCULATION 

Ii0.6,4HW/M2,3X,5HBETA=FI0.6) 

END!  It  lyOyOyOyOyly  lyOyOyOyOy  OyOyO) 
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MOSS 


CALCULATION  OF  t*  =  CONSTANT 


«  X  G  (,' 

OIMLNS ION  XLAMX ( 51 ) , T AUMXM ( 50 ) , AlPHA ( 50 , 50 ) , ZPRME ( 50 ,50 ) , FMAXM { 5 1 ) 
I » Y  L  A  M  X ( 5 II » TAUX! 51 ) ,ALPA( 51 ) ,ZPRMX( 50,50) 

ALFAF!X,Y)=(<<5./6. ) +(LOGF ( X )♦ { 5./3. ) »  LOGF  C Y) ) ) *SQRTF ( 1 . -Y** { -2  - /3 
1.  )  ) )/<  SORT F (-(LOGF {X ) )- { 5. / 3. ) *LOGF { Y ) ) ) 

T  A  U  M  X  F ( X ) =  (X*»{ 2./ 3. )+2. ) •SQRTF { X»* { 2./3. )-l. ) 

ZPRMF ( X,Y,Z  >=X*Z-( X*Y-YeSQRTF{ X**2.+! 10. /3. )•{ 1 ,-Y* • (~2. /3 •))))/( 2 
1.»S0IUF<  l.-Y»*(-2./3.m 

11  READ  INPUT  TAPE  5 , 100 , F MAX , XL AM , DE LFMX , DEL AM , M , N  . 
rMAAM(l)=PMAX 
XLAMXI 1 ) =XLAM 
DC  ?0  l  =  l ,  M 

WRITE  OUT  PUT  TAPE  6,  1 0 1 ,  I  ,  FMAXMI  I)  ,  DELFMK ,  DEL  AM 
DU  <0  J=i,N 

T  AUMXM ! J )  =  T  AUMXF ( XLAMX ( J ) ) 

ALPHA!  I, J  )=ALFAF< FMAXM!  I ), XLAMX! J)  ) 

ZPRME!  I , J )=ZPRMF( ALPHA!  I , J  )  , XL AMX ( J ) , T AUMXM ( J ) ) 

XLAMX < J+l )=XLAMX( J  M-DEL AM 

30  cuntinuf: 

WRITE  OUTPUT  TAPE  6, 1 02 , ( J , ALPHA ( I , J ) , ZPRME ! I , J ) , TAUMXM {  J ) , XL AMX ( J 

1  )  ,  J  =  1 ,  N  ) 

FMAXM!  1  +  1  )=FMAXM( I  I+DELFMX 
20  CONT  I NUE 

READ  INPUT  TAPE  5 ,  100 , X ALP A , YLAM ,DELALP , DE YLAM , L , K 
ALPA!  1  )=XALPA 
YLAMX (  1  )  =YLAM 

DO  50  l  =  1  *  L 

TAUX!  I  )=TAUMXF! YLAMX! I ) ) 

WRITE  OUTPUT  TAPE  6,203 

WRITE  OUTPUT  TAPE  6, 201 , I , TAUX (I) , YLAMX (I) , DELALP.DEYLAM 

DO  60  J  =  1,K 

ZPRMX (  I, J  )  = ZPRMF! ALPA! J ), YLAMX! I ) ,TAUX{ 1 ) ) 

ALP  A!  .3  +  1  )  =  ALPA  (  J  l+DELALP 
60  CONTINUE 

WRITE  UUTPUT  TAPE  6 , 202 , ( J , ALP A I J ) , ZPRMX { I , J ) , J  =  l , K ) 

,  YLAMX! I+l)=YLAMX( I I+DEYLAM 
50  CONTINUE 
GO  TO  11 

100  FORMAT  !4F 10,6, 216) 

101  FORMAT  ( 1H1,5X,2HI=I3,5X,6HFMAXM=F7.3,5X,7HDELFMX=F7.3,5X,6HDELAM= 
1F7. 3// ! 2X,  IHJ, 5X, 5HALPHA ,9X, 5HZPRME , 9X , 6HT A UMXM , 8X , 5HL AM AX ) ) 

102  FORMAT!  IX, I  2 , 3X , F 9 . 5 , 5X , F 1 0. 5 » 5X , F9. 4, 5X , F9 .4 ) 

201  FORMAT  ! 5X , 2H I  =  I  3 , 5X , 7HTAUMAX=F 7. 3 , 5X, 6HLAMAX=F 7^3 , 5 X , 9HDEL ALP 

1HA=F7.3,5X,6HDELAM=F7.3//! 2X , IHJ , 5X , 5HALPHA ,9X» 5HZPRME ) ) 

202  FORMAT  I  IX  ,  1 2 , 3X , F9 . 5 , 5X , F 10. 5 ) 

203  FORMAT! 1HU25X, 30HCALCULATI0N  OF  TAUMAX=CONST ANT/ { 26X, 29HAS  A  FUNC 
IT  I  ON  OF  ALPHA, ZPRIME ) ) 

END! 1 ,1,0,0, 0*0, 1, 1,0, 0,0, 0,0, 0,0) 
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T*  =  CONSTANT  WITH  CHECK  AND  TEST 


01  MENS  LON  ALPA( 51  )# YL AMXi 51) tTAUXI 51 ),ZPRME( 50, 50), CHECK  151,511 ,FM 
1AX!51,51),M151,51) 

ZPRMF!X,Y,Z5=X*Z-!  X» Y-Y*SQRTF  IX**2.+(10./3«  )*  (1  .-Y*»I-2./3.  ) )  )  )/  (2 
l.»S<)RTF!  U-Y**(-2./3.) ) ) 

TAU.tXF(X)-(X**«2./3.  ) +2. ) *SQRTF t  X»* ( 2. /3.  )-l.  ) 
FCTF(W,X,Y,Z)»!W»*(-5./3. )  MEXPF  <-! ( Z-X* Y) /W) **2. ) 

11  READ  [IMPUr  TAPE  5 , 100, XALPA , YLAM ,DELALP , DEYLAM, L , K 
Y=EXPF(-5./6. ) 

ALPa( l ) =XALPA 
YLAMX ( 1 ) =YLAM 
DO  50  1  =  l »  L 

T  AUX  ( 1 )=TAUMXF( YLAMX( I ) ) 

WRITE  OUTPUT  TAPE  6,203 

WRITE  UUTPUT  TAPE  6, 201, I ,TAUX(I) ,YLAMX(I) , DELALP , DEYLAM 

DU  60  -J  =  1,K 

ZPRME! I,  J)=ZPRMF!  ALPA ( J ) , YL AHX U  ),TAUX(I  )) 

FMAXl I,J )=FCTF ( YLAMX ( I ) , ALPA ( J ) , TAUX ( I ) , ZPRME ( I , J ) ) 

CHECK  ( I ,  J  )  3  FMAX(  I  ,J  )*(YLAMX( I )»*v5./3.  )  ) 

I F (CHECK (  1 ,  J  )- 1  •  )  51,51,52 

51  I F ( CHECK ( I , J ) -Y )  52,53,53 

52  M<I,J)=0 
GO  TU  60 

53  M  (  I  •  J  )  =  l 

60  ALPA<  J+l)=ALPA<  J  KDELALP 

WRITE  UUTPUT  TAPE  6,202,  < J , ALPA ( J ) , ZPRME ( I , J ) , FMAX I  I , J ) .CHECK (I , J 
l),M(I,J),J=l,K) 

50  YLAMXI  I  +  l  )=YLAMX(  IdDEYLAM 
GO  IU  11 

100  FORM AT (4F10. 6,216) 

203  FORMAT! IH1,25X,30HCALCULATIQN  OF  T A UMAX -CON ST ANT/ I26X, 29HAS  A  FUNC 
1 TION  OF  ALPHA, ZPRIME  )  ) 

201  FORMAT! 5X,2HI=I 3 , 5X , 7HTAUMAX=F 7. 3 , 5X , 6HLAMAX*F7. 3 , 5X ,9HDELALPHA=F7 
i.3,5X,6HDELAM*F7.3//(2X, 1H J , 5X , 5HAL PHA, 10X, 5HZPRME , 8X, 4HFMAX, 8X , 5H 
2CHECK.8X, 4HTEST )  ) 

202  FORMAT! IX, 1 2, 3X , F9.5 , 5X, F 10. 5 , 5X ,F9. 5 , 3X»F9.5 ,6X , II ) 

END! I, 1,0, 0,0,0, 1, 1, 0,0, 0,0, 0,0,0) 
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CALCULATION  OF  At 


DIMENSION  Z PRME ( 5 , 5 ) , ALPHA { 5 ) , XL AMX ( 5 ) , XLAMI 4 5 , 5 ) ,XL AM2 ( 5 , 5 ) , TAUX ( 
15),  T1!5,5),T2!6,5) ,FMAX(5,5  >,WIDTH(5,5) ,RESULTl5) ,FX2(5,5) ,FX1(5,5 
2) 

COMMON  ZPRME, ALPHA, F12, I, J 

FXXF!W,X,Y,Z)=vX*»(-5./3, )) «EXPF(-( { ( Z-W»Y ) /X ) **2 .  )  ) 

ZPKMF (X, Y,Z )=X*Z-(X*V-Y*SQRTF (X**2.*(10./3. )•( 1 .-Y«*!-2./3. ) )  )  )/(2 
l.*SQRTF( l.-Y»»(-2./3. ) )) 

TAUF( X)=IX**I2./3. )*2. ) *SQR TF { X*» C 2 . /3. 1-1. I 
2!)  READ  INPUI  f  APE  5, 602, EPS 

READ  INPUT  TAPE  5, 600, XALFA  ,XLAM,DE ALF ,DEXLM, M *N 
ALPHA ( 1 ) =X ALF A 
XLAMX I  1 ) =XL AM 
DO  90  I  =  L , M 

WRITE  OUTPUT  TAPE  6 , 601 , I , XLAMX ! I ) 

TAUX ( I)=TAUF(XLAMX(l) ) 

DO  91  J  =  l  ,  N 

ZPRME ( I  ,  J  )  -  ZPRMF(ALPHA(J),XLAMX(I ) , TAUX (I  )) 

FMAX(I,J)  =  FXXF(ALPHA{ J ) , XLAMX 1 1 ) , TAUX (I) , ZPRME ( I , J ) ) 
F12=0.5*FMAXt If J) 

DUMMY 

XLl=l. 00001 
XLK=XL AMX ( I ) 

XL2  =!5.0)*XLAMX{ I ) 

CALL  BRUUI ( XLM, XL 2, RESULT, EPS, DUMMY) 

IFIRESULT I  2 ) -  XL  M  )  50,51,50 
51  IF(RESULT(3)-XL2)  50,52,50 
52  WRITE  OUTPUT  TAPE  6,604 
GO  TO  91 

50  XLAM2I I »J)=RESJLT( 1  ) 

CALL  BROOTCXLI, XLM, RESULT, EPS, DUMMY) 

IFIRE SULT (2J-XL1)  40,41,40 
41  IFIRE SULT I 3J-XLM)  40,42,40 
42  WRITE  OUTPUT  1APE  6,604 
GO  TO  91 

40  XLAMI I  I , J )=RESULT( 1 ) 

TLI  I, J)=TAUF(XLAM1( I,J)  ) 

T2I  I, J)=TAUF{ XLAM2I I,J)  ) 

WIDTH  I  I , J ) =  T2I I , J)-Tl( 1, J) 

CHECK=  FMAX 1 1 , J ) *  I XLAMX 1 1 )**(5./3. ) ) 

Y=EXPFI-5./6. ) 

IF  (Y-CHECK)  20,20,21 
21  WRITE  OUTPUT  TAPE  6, 603, CHECK 
GO  TO  91 

20  IF  I  CHECK- 1 . )  30,  30,21 

30  FXU I , J ) =F  XXF ! ALPHA ( J ) , XL AM  1(1 ,J) , Tl ( I , J ) , ZPRME 1 1 , J ) ) 

FX2I  I , J )=FXXF< ALPHA! J),XLAM2( I , J) , T2 1 1 , J ) , ZPRME 1 1 , J ) ) 

WRITE  OUTPUT  TAPE  6 ,605 , J , ALPHA  I J ) , ZPRME (I , J ) , WIDTH U „ J ) , FX1 ( I , J ) , 
1FX2 ( I , J ) »  XLAMK  I , J ) , XL A M2  1 1  ,J),FMAX(I,J) 

91  ALPHA!  J  +  l  )=ALPHA(  JM-DEALF 
90  XLAMX l l+l )=XLAMX(l)*DEXLM 
GO  TO  25 

600  FORMAT  (4F10. 6,215) 

601  FORMAT! 1H1,5X»2HI=I2,5X, 6HL AMAX=F1 0. 6, / ,/(2X,lHJ,5X, 5H ALPHA, 10X , 5H 
1ZPRME , 9X, 5HWIDTH, 10X, 2HF 1 , 1  OX ,2HF2 , 9Xi 5HXLAM1 »9X»  5HXLAM2 , 9X , 4HFMAX 
2)) 

602  FORMAT  (F10.6) 

603  FORMAT! IX, 46HCHECK  DOES  NOT  MEET  PHYSICAL  CRITERION  CHECK=F10.6) 
604  FORMAT  I 31HN0  ROUT  OREVEN  NUMBER  OF  ROOTS) 

605  FORMAT! IX , I  2 , 4X , F 10. 6, 5X, FI 0. 6,4X , F 10. 6, 5X , F9. 6 ,3X, F9.6 ,3X * F 10. 6, 4 
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BISECTIONAL  ROOT  FINDER  FOR  At  PROGRAM 


SUBROUTINE  BR03T  ( X 1 , X2 *RESULT, C PS , DUMMY ) 

DIMENSION  ZPRMCIb, 5) .ALPHA! 5) .XLAMX 15) , XLAMl ( 5 , 5 ) ,X LAH2 ( 5, 5 ) , T AUX ( 
m,Tl(5«5),  T2(b*5)  ,  FMAX<5»5  )# WIDTH ! 5*5 ) .RESULT ( 5)  ,FX2  ( 5  *  5 )  *  FX1  <  5.  5 
2) 

TX1  =  X  1 
TX2=X2 

F1=UUMMY( TXl ) 

F2= DUMMY! 1X2) 

0  =  2 •»  EPS 
IF! f l ) 1,2, 2 

1  IF!F2)9,3,3 

2  IF(F2)3»9»9 

3  U=.5» ( TX I+TX2 ) 

F3=DUMMY ( U ) 

IF!F3)4,5,5 

4  I  F  (  F  I  )  7 , 6  »  6 

5  IF(FL)6,7,7 

6  F2=F3 
TX2  =  U 
GO  TO  8 

7  F 1  =  F  3 
TX  I  =U 

8  IF(ABSF(TX2-TX1)-Q)9,3,3 

9  RESULT! 1 ) = . 5* < T X 1+ T X2 ) 

RESULT! 2)=TX1 
RESULT! 3)=TX2 
RESULT(4)=F1 

RESULT! 5)=F2 

RETURN 

END 
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FUNCTION  STATEMENT  FOR  At  PROGRAM 


FUNCTION  DUMMY ( X  ) 

Cl. YEN SION  ZPRME(5f 5) , ALPHA (5 ), XLAMX(5) ,XLAMl(5i 5)  *XLAM2(5t5) ,TAUX( 
15) ,T1<  5,5) ,T2(5,5) ,FMAX(5,5) .WIDTH (5, 5) , RESULT ( 5 ) , F X2 ( 5 , 5 )  ,FX1<5*5 
2) 

COMMON  ZPRME,ALPHA,Fi2, 1, J 

DUMMY  =  F12*( X»*( 5./3.  )  )-EXPF(-( ( ( ZPRME ( I ♦ J ) -ALPHA ( J ) • ( ( X** ( 2 ./3 •  )  )■» 
12. ) *SCRTF ( ((X**(2./3.))-l.)) ) / X ) *• 2 • ) ) 

RETURN 

END (  ly  1, 0,0,0, Of  If  1,0, 0,0,0, 0,0,0) 


< 
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RECONSTRUCTION  OF  VOLTAGE  TRACES 


DIMENSION  AC 50) ,Y( 50 ) ,XLAX(51 ) , F ( 50 ) , T AUX ( 5 0 ) »V(50) ,  T ( 50 ) , DENS ( 50 ) 
1 ,V£L<  5C  ) 

FCTF(W, X,Y,Z )  =  t W**  C-5./3. ) )*EXPF  t-( CZ-X*Y)/W)**2.  ) 
TAUF(X)=((X»*(2./3.))+2. ) *  SQR  T  F  C  ( X** ( 2. /3. ) )-l .  )  , 

92  READ  INPUT  TAPE  5, 900, ALPHA, XLAM ,DX, XLAMO, ZO, N 
READ  INPUT  TAPE  5 , 90 1 , RO , R 1 , R 3  ,  B  0 »  RC 
READ  INPUT  TAPE  5 , 902 , VO , XM , XN , X  I 
PI  =  3. 1415927 
XMu  =  P  I  »4.  et-=07 

C=  ( XMU*XN*XM*V0)/(6.0*(PI ** ( 3./2. ) ) *ALPHA*XI *( R0**2.  ) *(B0»*2. ) ) 

91  READ  INPUT  TAPE  5 , 903 , Z , ZLAST 
ZPRME  = ( Z-ZO ) / ( XLAMO ) 

XLAX( 1 ) =  X  L  A  M 

DO  90  1  =  1,  N 

T  AUX (  I  )  -T  AUF ( XL  AX (  I  )  ) 

F( I  )=FCTF(XLAX( I  )  t ALPHA*  T  «UX ( I  ), ZPRME) 

A(I ) =C  *F (  I  ) 

Y  (  I  '  =  SQRTF(  (A(  I )*«2. )  +  l • ) —A (  I  ) 

VI  I  )=  (  (P  I*BO)/RC )*( l.-Y( I ) )*( (R0**2. )/Y( I )  )*(  (  (R3**2. )-(ftl**2. ) )/ 
1 ( ( R3  *  *  2 • )-(  ( R0*»2 •  )  / Y(  I ) ) )  ) 

D  E  N  S (  I ) =  (  (XN/ (  (PI**(3./2. ) )* XLAMO* (R0**2. )  ) ) *Y( I ) *EXPF (-( ( ZPRME-A 
1LPHA*T  AUX  (  ())/XLAX(m*»2.)  }/XLAX(  I  ) 

VEL  (  I )  =  (  XIAMO/X I ) » (ALPHA* ( ( ZPRME-ALPHA* TAUX ( I) ) /XLAX ( I ) )*SORTF( 1.- 
1 ( XL AX ( I ) *  * ( -2 . / 3  .  )  )  )  ) 

T ( I ) =X I *TAUX ( I  ) 

90  X L A X (  1  +  1 )  =  XLAX(  I  ) +DX 

WRITE  OUTPUT  TAPE  6  *  904 , Z , ALPHA , XM f XN , ZO , VO , XLAMO , BO  ,X I 
WRITE  OUTPUT  TAPE  6,906 

WRITE  OUTPUT  TAPE  6, 905,  (  I  ,  V  (  I  )  *  T  (  I  )  ,  A  (  I  )  ,  T  AUX  CD  t  XL  AX  (  I  )  ,  DENS  (  I  )  , 
IVEL (  I  ) , 1  =  1, N  ) 

IF(Z-ZLAST)  91,92,92 

900  FORMAT  (5FL0.6,I5) 

901  FORMAT  (4F10.6,E10.4) 

902  FORMAT  (4E10.4) 

903  FORMAT  (2F10.6) 

904  FORMAT  (  IH 1 , 5X ,  2HZ  =  F5 . 3 , 7H  ME T.E RS  ,  3 X ,  6H AL PHA=F6.  3 , 3X  ,  5HMASS  =  E  10 . 3 , 6 
1H  KGRMS,3X,5HIONS=E10.3/(6X,3HZO=F5.3,7H  METERS » 3X , 3HV0=E10 . 2 , 6H  M 
2/ SEC, 3X,6HXLAM0=F5.3,7H  METERS ,3X,3HB0=F 7.4, 5H  W/M2 ,3X ,3HXI=E10. 3, 
34H  SEC,//  )  ) 

905  FORMAT ( 2X , 1 2 , 3X , F6 . 3 , 3X , E 10 . 3 , 4X , F7. 3 , 4X , F7 . 3 , 4X , F7 . 4 , 3X , E 10. 3 , 4X  , 
1E10.3) 

906  FORMAT (3X, 1H I , 6X , 1 HV , 1 IX , 1HT , 1 l X , lHA , 9X , 3HT AU , 5X , 6H L AMBDA ,7X , 4HDEN 
IS,9X,3HV'EL  ) 

END ( 1, 1,0, 0,0, 0,1, 1,0, 0,0, 0,0, 0,0) 
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